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The transport phenomena of Lennard- Jones molecules through a structural bottleneck 
driven by an external force are investigated by molecular dynamics simulations. We observe 
two distinct molecular flow regimes distinguished by a critical external force F c and find scal- 
ing behaviors between external forces and flow rates. Below the threshold F c , molecules are 
essentially stuck in the bottleneck due to the attractive interaction between the molecules, 
while above F c , molecules can smoothly move in the pipe. A critical flow rate q c correspond- 
ing to F c satisfies a simple relationship with angles, and the value of q c can be estimated 
from a simple argument. We further clarify the role of the temperature dependence in the 
molecular flows through the bottleneck. 

§1. Introduction 

Transport phenomena through a structural bottleneck have been widely observed 
and have attracted attention owing to the development of nanoscience. Many in- 
teresting bottleneck features have been observed and studied in granular flowsP® 
traffic flows,®®'^ evacuation processes,®'^® water flow through porous grains,^ 
microfluidicsp^ Laval nozzles, 11 -' and the Internet.^ Recent technical advances 
have enabled us to study the nanoscale transport system such as ion transport 
through protein channels in a lipid membrane artificial nanoporesp^ car- 
bon nanotubesp^ Knudsen gasesp^J'^ and domain wall motion in a wire with a 
notch ,[2TJ,[20j) r j i ] luS) a comprehensive understanding of transport phenomena 
through a bottleneck is essential. 

Unlike bulk systems, the bottleneck alters the transport properties. Jamming 
phenomena relating the formation of a permanent arch at the bottleneck^ and 
intermittent outflows in a funnel with various angles^ have been observed and in- 
tensively studied!^* In electrical transport through quantum wires, the effect of 
the geometry of the wire on the transmission coefficient and conductance has been 
an important issueP^'^ Recently, experiment and simulation studies have been 
performed to understand the structure and dynamics of water in a confined carbon 
nanotube, where it is discussed that confinements induce the phase transition not ob- 
served in a bulk systemP^*^ In the transport phenomena through the bottleneck, 
these collimator effects such as the curvatures and cross sections of the passage play 
a crucial role in the transport phenomena. In such situations, a useful approach to 
treat transport phenomena through a bottleneck is computer simulations. Various 
simulation techniques, such as driven lattice gasesj®'® cellular automata,^ and 
molecular dynamics (MD) 34) have been designed and applied to probe these trans- 
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port mechanisms. Numerically, the existence of scaling behavior of escape dynamics 
through a bottleneck^ ® and the power spectra of the flow in an hourglass 37 ' have 
been investigated. 

In this paper, we address the problem of the steady particle flow in a pinched 
pipe under an external force. We restrict our attention to a highly simplified two- 
dimensional model, considering geometric effects in the hope to gain some generic 
physical insight. We intend to evaluate the importance of the curvature of the pipe, 
particle area fractions, and external forces. This goal may be achieved through 
detailed simulations of the motion of particles. For these reasons, we perform MD 
simulations to investigate molecular flows in the pipe with various angles. 

The present article is organized as follows. In £[2] we discuss the model and MD 
simulation techniques. In fj3j the results of the MD simulation and data analyses 
are presented. Special emphasis is placed on the relationship between stationary 
molecular flows and external forces at various angles and scaling behavior of these 
relations. Then, we demonstrate that flow rates against angles give a simple relation 
and estimate the value of the critical flow rate. We investigate the fluctuation of 
the flow and the fundamental diagram for molecular flows in analogy to a traffic 
model. We further investigate the temperature dependence of the molecular flow 
and how flow rates depend on the ratio between molecular size and pipe width. For 
comparison, we observe a molecular flow that has only repulsive interactions. We 
conclude and discuss our results in 21 

§2. Model 

Let us consider a two-dimensional pinched pipe shown in Fig. []Ja). The pipe 
has symmetry configuration, like an hour glass, and is represented by the hyperbolic 
tangent y = 36 tanh(x tan 9°/36) with the angle 9 defined by the slope of hyperbolic 
tangent. 




(a) (b) 

Fig. 1. (a) Schematic of a two-dimensional symmetric pipe with angle 8, which is characterized 
by the slope of y — 36tanh(a;tan#/36). The rr-directional pipe length is Lx — 300<r and 
the narrowest pipe width is w. The wall is constructed of particles that are attached to other 
particles by a linear spring, and the y-directional boundary is the thermal wall with temperature 
T = O.le. (b) Schematic of left region of the pipe with angle 9 = 50°, 70°, 80°, 90°. 
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The molecules are interacting via Lennard-Jones (LJ) potential: 

'-(£>}■ (21) 

where e is the energy scale of the interaction, a is the diameter of each molecule, and 
r{j is the separation between particles i and j. The pair interactions are truncated 
at r^ = 3a for numerical convergences. The linear x-directional pipe length is 
L x = 300ct. In most of our simulations, the narrowest pipe width w is 3a, but we 
will later investigate the width dependences on the flow rate. The wall is composed of 
molecules that are attached to other molecules by linear springs with spring constant 
k = 20000e/o" 2 , but the molecules at both edges are fixed. The wall molecules are 
interacting with driven molecules via the same LJ potential. In the presence of the 
x-directional external force F, the current of the molecules will be maintained, and 
hence, the system will always remain in a nonequilibrium stationary state. 

In the horizontal boundaries, we adopt the periodic boundary conditions for the 
positions of molecules and the thermal walls for the velocity of molecules, which give 
the random velocity whose probability distribution function is given by 

1 ( 771 \ ^ mv^ 

m = ^m\kSf) Me ~' kBr (2 ' 2) 

with temperature T = O.le. Owing to the presence of the x-directional external 
forces, molecules keep flowing from left to right through the bottleneck. 

For comparisons of LJ molecular flows, we adopt the purely repulsive Weeks- 
Chandler- Andersen (WCA) potential given by 



0WCA(ry) = { ~ ' 4 J ' V (2-3) 

, otherwise 

which we call WCA for later discussion. 

MD simulation is carried out for the dynamics of the molecules under the differ- 
ent geometrical pipes. The equation of motion is solved by using the velocity- Ver let 
method with a time step At = O.OOIt, where r = a(m/ e) 1 / 2 . The mass of the 
driven molecule is m and that of the wall is 10m to avoid the change in the pipe 
configuration due to the molecular flow through the pipe. The initial positions of 
molecules are arranged not to contact each other, and the initial velocities of the 
molecules are given by the Maxwell distribution with temperature T = O.le. In our 
reduced units, the unit length a, unit energy e, molecular mass m, and Boltzmann 
constant /cb are set equal to unity. As shown in Fig. QJb), we vary the pipe angle 
= 50°, 55°, 60°, 65°, 70°, 75°, 80°, 85°, 90°. As the angles become large, each pipe 
represents a narrow straight passage with a narrowed entrance and constant pipe 
width, whereas the pipe with a small angle exhibits a pinched configuration. The 
pipe configuration less than angle 6 = 50° is not used since the y-directional pipe 
length becomes markedly short in our model. 
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To investigate molecular flows in the pinched pipe, we mainly measure the flow 
rate q in the pipe with various angles, which is the numbers of crossing-driven 
molecules during the time at the middle cross section. In the steady flow, the flow 
rate q should be a constant in any cross section owing to the continuity equation in 
the stationary state. After t = 1.5 x 10 4 r, we collect the statistical quantities as the 
system reaches a steady state. 

§3. Steady molecular flow 

3.1. Steady molecular flow at various angles 
30 
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Fig. 2. Plots of the flow rates q against the various external forces F in the pipe with angle 9 = 50°, 
55°, 60°, 70°, 80°, 85°, 90°. 

We study steady molecular flows through the pipes with various angles. The 
narrowest width of each pipe is three times the molecule diameter. The driven 
molecule number in the pipe is N = 3510. Figure [2] plots the flow rate q against the 
external force F under various angles. In the case of small external forces, driven 
molecules jam in the pipe and stick on the wall since the molecular interactions 
are sufficiently strong compared with external forces, while as the external forces 
increase, the nonlinear flow growth is observed and jamming is insignificant owing 
to the dominance of the external forces. These transport phenomena are clearly 
observed in the pipes with angle below 9 = 80°. On the other hand, as Fig. [3] shows, 
the distinct relation between the flow rate and external force is difficult to observe 
at the angle above 9 = 80° owing to the small amount of molecular flow caused 
by the markedly narrow pipe entrance and the passage at these angles. At angle 
6 = 90°, the molecular arch is observed in the low external force. Figure 0] exhibits 
the existence of a molecular arch in the entrance of the pipe with angle 9 = 90° 
under the external force F = 0.010e/o". This phenomenon could be caused by the 
straight corner at the pipe entrance. 

To investigate the molecular flow properties further, we perform the semilog plot 
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Fig. 3. Plots show the flow rates q against the external forces F at the angle 6 = 80° and 85° with 
variance. 



Fig. 4. A schematic pipe graph at the angle 6 = 90° and the external force F = O.OlOe/a. The 
entrance of the pipe exhibits arch formation. 



of the flow rate, log q, against the external force F. The inset of Fig. [5] indicates 
the existence of two distinct molecular flow regimes in the pipe with angle 6 = 60°. 
The flow rate and external force at the crossover of these regimes, q c and F c , are 
considered to correspond to the balance of external forces and molecular interactions. 
These transport phenomena are observed at the angle below 6 = 80° and q c and F c 
are determined at each angle. At the angle above 9 = 80°, it is difficult to distinguish 
characteristic crossover from the relation between flow rates and external forces. By 
using q c and F c in each pipe, the scaled flow q/q c against the scaled external force 
F/F c can be plotted in Fig. The solid and dotted lines are fitted by curves 
q/q c = exp(aF/F c — b) with a,b = (8.45,8.48) for molecular flow below F c and 
(0.445, 0.431) for molecular flow above F c , respectively. Thus, it is shown that flow 
rates as a function of external forces are given by the two distinct simple exponential 
forms and are independent of the pipe geometry in these regimes. 

The angle dependences on q c can be simply studied. Figure [6] is the plot of 
q c as a function of the angle. At a large angle, pipe configurations represent the 
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Fig. 5. Plot of the scaled flow q/q c against the scaled external force F/F c at the angle 8 = 50°, 
55°, 60°, 70°, 80°. The solid line is fitting curve q/q c = exp(aF/F c - b) with a, b = (8.45, 8.48) 
for molecular flow below F c and the dotted line is a, b = (0.445, 0.431) for molecular flow above 
F c . The inset shows a semilog plot of the flow rate and the external force in the pipe with angle 
9 = 60°. We define the crossover of two flow lines as F c and q c . 
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Fig. 6. Critical flow rate q c as a function of the pipe angle 9. The solid line is the fitted curve 
a cos 2 9 + qo with a = 25.9(e/m) 1//2 /(T and go = 4.21(e/m) 1//2 /cr. The unit of the vertical axis is 
(e/m) 1 / 2 /a. The inset shows the plot of the critical external force F c [e/a] against the angle 0. 

reduced passage so that it is plausible for the molecular flow to decrease against the 
angles. We are aiming at a rough theoretical estimation of variations of q c against 
the angles. In the presence of F c , the system is considered to be balanced between 
the external force and the molecular interaction. We assume that the practical 
external force along the pipe with angle 6, F c cos6, drives molecules in the pipe. 
This practical external force is set off in the y-direction owing to the symmetric pipe 
configuration; thus, the x-component of the practical external force F ~ F c cos 2 9 
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could be dominant for molecular flows. Since we consider the steady state, the 
particle density in the filled pipe region p and the terminal velocity v ~ \xF with the 
mobility \i are considered to be constant. Thus, the flow rate at the balanced point 
is approximated by 

q c ~ a cos 2 9 + qo, (3T) 

with fitting parameters a and qo, where qo is the extrapolation of the critical flow 
rate at 9 = 90°. In Fig. El the solid line represents (13-li with a = 25.9(e/m) 1 / 2 /cr 
and qo = 4.21(e/m) 1 / 2 /cr and is well fitted at the angle below 9 = 80°. On the 
other hand, the relation F c and the angle are almost independent of the angle under 
9 < 75° in the inset of Fig. [6) At the angle above 9 = 80°, the value of F c that 
corresponds to the value of q c cannot be determined since the crossover is not clearly 
observed. 

We further estimate the value of q c , which corresponds to F c . The critical flow 
rate q c is given by 

q c = nv c S, (3-2) 

where n = 4</>/(7r<7 2 ) is the number density with the area fraction cj), v c is the critical 
velocity, and S is the cross section. The critical velocity is represented by the escape 
velocity from the bottom of LJ potential as v c = ^/2\(pLj(a)\/m with the distance 
at the minimum potential a = 2 1 / 6 cr. Thus, the critical flow rate can be estimated 
by 

q c = 7=<P- ( 3-3 J 

7r yjma 

The area fraction at the region where the particles accumulate is eft ~ 0.7 in our 
model, so that we can estimate the value of the critical flow rate as q c ~ 3.8^ e/m/a. 
This value is close to the critical flow rate qo at the angle 9 = 90°. Thus, these 
theoretical estimations are plausible for the explanation of our molecular transport. 

3.2. Fluctuation of the flow 

We next investigate how the fluctuation of the flow depends on the external 
force. To clarify the characteristics of the fluctuation, we introduce the normalized 
standard deviation a* obtained by the linear interpolation of the standard deviations 
cr[qi/q c ] and crfa/qc], where q\ and q<i are the closest points from q c (qi < q c < q-i). 
Figure [7] shows the standard deviation of the scaled flow rate cr[q/q c ]/a* against 
the scaled external force F/F c at the angle 9 = 50°, 55°, 60°, 65°, 70°, 80°. As 
Fig. [7] shows, the scaled standard deviation a[q/q c ]/a* increases gradually below the 
external force F < F c , where the particles are essentially stuck in the bottleneck 
and the fluctuations are suppressed. For the large external forces F > F c the scaled 
standard deviation is almost constant in any pipe angle and the fluctuation of the flow 
is insensitive to the external force. These fluctuation phenomena almost correspond 
to Fig. [5] (the relation between the scaled flow rate and the scaled external force). It 
is interesting that the fluctuation of the flow does not have any specific peak around 
F c , and the critical fluctuation continues above F > F c , while we do not make clear 
its reason. 
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Furthermore, we investigate how the fluctuation depends on the angle. Figure. 
[8] shows the standard deviation of the flow a[q] against the angle 9 at the constant 
external forces F = 0.020, 0.016, 0.004e/cr, and the inset of Fig. [8] shows the cor- 
responding flow rate q against the angle 9. In each external force, the standard 
deviation of the flow rate is substantially constant and does not exhibit any fluctu- 
ation peak against the angle. Although we have investigated the relation between 
the fluctuation and the angle, the angle dependence of the fluctuation of the flow is 
not observed in our model. 
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Fig. 7. Standard deviation of the scaled flow rate a[q/q c ]/a* against the scaled external force F/F c 
at the angle 9 — 50°, 55°, 60°, 65°, 70°, 80°. The normalized standard deviation a* is the linear 
interpolation of the standard deviations a[qi/q c ] and a[q2/q c ], where qi and Q2 are the closest 
points from q c (qi < q c < 52) in each angle. 
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Fig. 8. Standard deviation of the flow rate a[q] against the angle 9 in the external force F = 
0.020, 0.016, 0.006e/ct. The inset shows the corresponding flow rate q against the angle 9. 
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3.3. Fundamental diagram 

For traffic flows, the functional relations between vehicle currents and vehicle 
densities have attracted much attention, which is called the fundamental diagram.® 
In our model, we investigate the relation between the molecular flow rate and the 
molecular area fraction under constant external forces. We define the area fraction (ft 
as the average area fraction in upstream regions (0 < x < 150cr), which corresponds 
to the accumulated molecular height from the middle of the pipe. Figure [9] shows 
the scaled flow rate q/q c as a function of (ft in the pipe with angle 9 = 60° and 
the narrowest pipe width 3a. In each scaled external forces F/F c , scaled flow rates 
exhibit the linear dependence of (ft i n l° w area fractions, while scaled flow rates 
decrease in high area fractions. These behaviors in our system are independent of 
F c and similar to traffic flow behaviors. 

The existence of the metastable branch is often discussed in the fundamental 
diagram of the traffic flow, while the present model does not exhibit such behaviors. 
As the previous subsection has shown, we cannot observe any fluctuation peak of 
the flow in both Figs. [7|and[8l These results likely relate the lack of the metastable 
branch in our model. Thus, the flow rate may be represented by a single curve in 
our fundamental diagram. 

3.4. Width dependence of flow rate 

In the situation of gravity-driven granular flow, jamming transitions depend on 
the ratio between the exit size and the grain sizeP^ Thus, we investigate the relation 
between flow rates and pipe widths under constant external forces. We fix a pipe 
configuration with the angle 9 = 60°, but vary the pipe width of the narrowest part 
from w = 1.5a to w = 6a. The driven particle number is N = 3510 and the external 
force is F = 0.012e/a, which is above F c and F = 0.008e/<7 below F c . Figure [TOl 
shows the flow rate q against the narrowest pipe width w. In F = 0.012e/o", no 
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Fig. 10. Flow rate q as a lunction of the pipe width at the angle 
F = 0.012e/<r and 0.008e/cr. Dotted lines are connected data. 
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molecular flow is observed in the pipe width under w = 2a, which is too narrow 
to pass through molecules. In F = 0.008e/<7, we observe no molecular flow below 
w = 4.5a owing to the low external forces. On the other hand, there exists the 
flow when the width is larger than w = 2a for F = 0.012e/<t and w = 4.5cx for 
F = 0.008e/<7. Although the proportional constant depends on the external force, 
we observe that q is almost proportional to the cross section of the pipe in this flow 
region. 

3.5. Temperature dependence 

To investigate the temperature dependence on the flow rate, we consider LJ 
molecule at the temperature T = 0.1, 0.2, 0.3, 0.4e in the pipe with the angle 9 = 60° 
and the width 3a and molecular number N = 3510. Figure [TT] shows the scaled flow 
rate q/qF, where qF is the flow rate at the external force F = 0.020e/<r, as a function 
of the external forces. As we have shown in Fig. [5] there is the "jammed region" 
for F < F c at the temperature T = O.le. However, as the temperature becomes 
large, LJ molecules can flow smoothly in the "jammed" region (F < F c ) and the 
flow rate increases linearly against the external force. It is known that T = 0.5e 
is the transition temperature from the gas phase to the liquid phase for the wide 
range of the density.^ As Fig. [11] shows, we can observe the smooth flow at the 
temperature T <; 0.4e, so that this "jamming" phenomena likely relate gas-liquid 
phase transition. At the low temperature T = O.le, LJ molecules are regarded as 
in liquid phase, so that the correlated motion of the liquid molecule may cause no 
molecular flow and dominate the "jamming" phenomena in the pipe. 

We also investigate the temperature dependence on the flow in the constant 
external forces. Figure [12] shows the flow rate q with the standard deviation as a 
function of temperature T in the external field F = 0.004e/<7 and 0.012e/<7. In 
the case of the external force F = 0.004e/o", LJ molecules are almost stuck at the 
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Fig. 11. Scaled flow rate q/qF as a function of the external force in the pipe with angle 6 
60°, where qp is the flow rate at F = 0.020e/<r, for LJ molecule at the temperature T 
0.1, 0.2, 0.3, 0.4e and WCA molecule at the temperature T = O.le. 



low temperature and begin to flow above about T = 0.2e. In the case of the large 
external force F = 0.012e/<r, LJ molecules flow even in the low temperature and the 
"jamming" phenomena do not happen. At high temperature, "jamming" phenomena 
do not occur and, in addition, we have observed that "jamming" phenomena depend 
on the external force. 
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Fig. 12. Flow rate q with the standard deviation as a function of the temperature T at the external 
field F = 0.004e/a and 0.012e/a in the pipe with angle 9 = 60°. 
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3.6. WCA potential 

Thus far, we have investigated the transport of LJ molecules that have both 
attractive and repulsive interactions against other molecules. To extract the role of 
attractive interactions, we study the molecular flow with WCA potential <^>wca in eq. 
(|2-3p . As the same condition of LJ molecular flow, the angle 9 = 60°, the width 3a, 
and molecular number N = 3510 are ued. The plot of the triangle in Fig. [TT1 shows 
WCA molecule flow. LJ molecular flow at the temperature T = O.le is represented 
by two exponential functions as mentioned earlier, while WCA molecular flow can 
be fitted using a single exponential function. In the large external forces, LJ and 
WCA molecular flows exhibit an identical behavior, since external force is dominant 
for both flows. 

Furthermore, we investigate whether the arch effect appears for the WCA model 
at the pipe with the angle 9 = 90°. The WCA molecule can flow smoothly in the 
pipe in any external force and even at low temperature, so that we observe no arch 
effect for the WCA model. These results indicate that the attractive force and low 
temperature are dominant for "jamming" phenomena. 

§4. Discussion and conclusions 

We have studied the molecular transport through bottlenecks, which is inspired 
by nanoscale transport phenomena such as biological ion transport. In this study, 
we simplify the bottleneck as the two-dimensional pinched pipe characterized by the 
angle 9 and width uu. As the external forces increase, we can observe two distinct 
molecular flow regimes. Using the value at the crossover of these regimes, which is 
considered to correspond to the balance of the molecular interaction and the external 
force, the scaling behaviors of the flow rate and external force are observed in our 
model. These scaling behaviors show that the flow rate does not depend on the 
pipe configurations at the angle below 9 = 80°. However, at the large angle, we 
cannot clarify the relation between the flow rate and external force owing to the 
small amount of molecular flow. The flow rate as a function of the angle is also 
observed, and this relation is verified using a simple theory. We further exhibit 
the value of the critical flow rate from rough theoretical estimations. These double 
exponential forms of the flow rate against the external force can be clearly observed 
at the low temperature T = O.le. On the other hand, as the temperature increases, 
"jamming" region disappears even in the small external force and this "jamming" 
phenomena are likely related to the gas-liquid transition of the 2D LJ phase diagram. 
In the fundamental diagram, the change in the fluctuation in the metastable region 
is often discussed, whereas we cannot observe any fluctuation peak in our model. 
These results likely relate to the lack of the metastable branch and thus the flow 
rate may be represented by a single curve in our fundamental diagram. We also 
investigated the WCA molecular flow to extract the role of the attractive forces and 
found that attractive forces at the low temperature are dominant for the "jamming" 
phenomena. 

The extension of this model to other transport models is straightforward. As 
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examples, grand canonical molecular simulations have been developed to examine 
the molecular transport by introducing explicit chemical potential gradients into the 
systemPJJH) 

As emphasized in the introduction, we study a very simple model of the pipe 
with a bottleneck, while our model qualitatively shows the important features of the 
molecular transport through the bottleneck and is of importance for the architecture 
and functionality in the emerging field of nanomaterial systems. 
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